paquetes
library (DT)
library (dplyr)
library (terra)
library (raster)
library (leaflet)
library (sf)
library (ggplot2)
library (plotly)
capas
provincias <-
st_read(
"https://github.com/tpb728O-programaciongeoespacialr/2021ii/raw/main/datos/ign/delimitacion-territorial-administrativa/provincias.geojson",
quiet = TRUE
)
provincias <-
provincias %>%
st_transform(4326)
cuencas <-
st_read(
"https://raw.githubusercontent.com/MaureenArg/datostarea/master/cuencas.geojson",
quiet = TRUE
)
altitud <-
rast(
"/vsicurl/https://raw.githubusercontent.com/tpb728O-programaciongeoespacialr/2021ii/master/datos/worldclim/altitud.tif"
)
denu_depu <-
st_read(
"/vsicurl/https://raw.githubusercontent.com/MaureenArg/datostarea/master/denu_depu.csv",
options = c(
"X_POSSIBLE_NAMES=decimalLongitude",
"Y_POSSIBLE_NAMES=decimalLatitude"
),
quiet = TRUE
)
st_crs (denu_depu) <- 4326
tabla
denu_depu %>%
st_drop_geometry()%>%
select (TIPO_den, año, TIPO_inf, CATEGORIA_) %>%
datatable (rownames= FALSE, filter= "top", class= "hover cell-border stripe", caption= "Tabla 1: Tipo de denuncias ambientales según año", colnames = c( "Tipo de denuncia", "Año", "Detalle de la denuncia", "Categoría"), extensions = "Buttons", options = list (language = list (url = "//cdn.datatables.net/plug-ins/1.10.11/i18n/Spanish.json"),dom = "Bfrtip", buttons = c("csv", "excel", "pdf")))
mapa
altitud_rl <- raster(altitud)
leaflet() %>%
addTiles (group = "OSM") %>%
addProviderTiles("Esri", group = "Esri") %>%
addProviderTiles("CartoDB", group = "CartoDB") %>%
addRasterImage(
altitud_rl,
opacity = 0.6
) %>%
addPolygons(
data = cuencas,group = "Cuencas hidrográficas",
color = "black",
fillColor = "blue",
stroke = TRUE,
weight = 1.0, popup = ~NOMBRE,
) %>%
addPolygons(
data = provincias,group = "Provincias",
color = "black",
fillColor = "Transparent",
stroke = TRUE, popup = ~provincia,
weight = 1.0
) %>%
addCircleMarkers(
data = denu_depu,
stroke = F,
radius = 4,
fillColor = 'red',
fillOpacity = 0.5, popup = ~TIPO_inf,
) %>%
addLayersControl(baseGroups = c ("OSM", "Esri", "CartoDB", "Cuencas hidrográficas", "Provincias"))
## Warning in validateCoords(lng, lat, funcName): Data contains 1 rows with either
## missing or invalid lat/lon values and will be ignored
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded ellps WGS 84 in Proj4 definition: +proj=merc +a=6378137
## +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=m +nadgrids=@null
## +wktext +no_defs +type=crs
## Warning in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj =
## prefer_proj): Discarded datum World Geodetic System 1984 in Proj4 definition